Tomography of high-energy nuclear colhsions with photon-hadron correlations 
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Within the next-to-leading order (NLO) perturbative QCD (pQCD) parton model, suppression 
of away-side hadron spectra associated with a high pr photon due to parton energy loss is studied 
in high-energy heavy-ion collisions. Dictated by the shape of the 7-associated jet spectrum in NLO 
pQCD, hadron spectra at large zt = Pt/Pt ~ ^ more sensitive to parton energy loss and 
therefore are dominated by surface emission of 7-associated jets, whereas small zt hadrons mainly 
come from fragmentation of jets with reduced energy which is controlled by the volume emission. 
These lead to different centrality dependence of the 7- hadron suppression for different values of zt- 
Therefore, a complete measurement of the suppression of 7-triggered hadron spectra, including its 
dependence on the orientation of the 7-hadron pair with respect to the reaction plane, allows the 
extraction of the spatial distribution of jet quenching parameters, achieving a true tomographic 
study of the quark-gluon plasma in high-energy heavy-ion collisions. 

PACS numbers: 12.38.Mh, 24.85. +p; 25.75.-q 



Jet quenching [T] or suppression of large pT hadrons, 
caused by parton energy loss due to strong interaction 
between jets and the dense medium, has become a pow- 
erful tool for the study of properties of the quark-gluon 
plasma [2 in high-energy heavy-ion collisions. Exper- 
imental studies of jet quenching include suppression of 
single [31 dihadron (back-to-back) |Hj and 7-hadron 
spectra. The strong suppression of large pt single hadron 
spectra and its centrality dependence at the Relativistic 
Heavy-ion CoUider (RHIC) [31 [5], indicate a picture of 
surface emission of jets. Most of jets, initially produced 
at the center of collisions, undergo large amount of energy 
loss and do not contribute to the final observed large pT 
hadron spectra. High px dihadrons, on the other hand, 
come not only from jet pairs close and tangential to the 
surface of the dense medium but also from punch-through 
jets originating from the center of the system [6], leading 
to their increased sensitivity to the initial gluon density 
as compared to single hadrons. 

We will focus in this paper on the study of 7- 
triggered away-side hadron spectra in heavy-ion colli- 
sions within the next-to-leading order (NLO) perturba- 
tive QCD (pQCD) parton model. Since a photon does 
not interact in QCD with the dense medium, its en- 
ergy approximately reflects that of the initial jet in 7-jet 
events prior to jet propagation through the medium. One 
can therefore study the medium modification of the full 
jet fragmentation function (FF) [3 |8]. By selecting 7- 
hadron pairs with different values of = PtIPti which 
could be larger than 1 due to radiative correction in NLO 
pQCD, one can effectively control hadron emission from 
different regions of the dense medium and therefore ex- 
tract the corresponding jet quenching parameters. 

Within pQCD parton model, the NLO corrections 
(aeO^s) to photon and photon-hadron production cross 
sections include 1-loop corrections to 2 ^ 2 tree level 



processes and 2^3 tree level contributions which have 
two-cutoff parameters, 5s and 5c, to isolate the soft 
and coUinear divergences in the squared matrix elements 
[HI Hn] . This results in a set of two-body and three-body 
weights that depend on the cut-offs. The dependence 
cancels when the weights are combined in the calcula- 
tion of physical observable, such as inclusive photon and 
photon-hadron cross sections. In this paper, we will use 
a Monte Carlo implementation [T^ of the NLO pQCD 
calculation for the invariant cross section of photon and 
photon-hadron production. 

Photon production in pQCD includes direct photons 
from the Compton and annihilation subprocesses and 
bremsstrahlung photons from jet fragmentation in high- 
energy p + p collisions. Since the bremsstrahlung pho- 
tons are always accompanied by nearly coUinear hadrons, 
an isolation-cut can be applied on the photon signal to 
suppress bremsstrahlung-like photons both in theory [TD] 
and experiments Even though such isolation-cut 

cannot eliminate fragmentation contributions, it can sig- 
nificantly reduce the fraction of fragmentation photons 
to about less than 10% for p^ > 7 GeV/c. For the 
study of 7-hadron correlation in this paper, we will fo- 
cus mainly on photons with isolation cuts. Therefore, we 
can neglect those photons that are produced via induced 
bremsstrahlung |12j . jet-photon conversion [13 and ther- 
mal production |14l I15j in high-energy heavy-ion colli- 
sions. 

To take into account parton energy loss in 7-hadron 
correlation in high-energy heavy-ion collisions, we will 
use an effective form of medium modified FF's as in pre- 
vious studies [6l [16] with an average energy loss, 
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for a quark (gluon energy loss is 9/4 of a quark) produced 
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at a position r and traveling along the direction n, where, 
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dL 

is the parameterized |17| average quark energy loss per 
unit length in a 1-d expanding medium with an initial 
average gluon density pq at time tq. The parameter eo 
is the initial quark energy loss per unit distance and is 
proportional to pq. A hard-sphere nuclear overlap geom- 
etry is used as in previous studies [B] [16], which differs 
at most about 10% from a Wood-Saxon geometry. 

We will use the KKP parameterization (TH] for parton 
FF's in vacuum. The parton distributions per nucleon in- 
side a nucleus are assumed to be factorizable into parton 
distributions in a nucleon given by the CTEQ6M param- 
eterization fJW and a nuclear modification factor [20], 
including the isospin dependence. In the following, we 
will use values of the energy loss parameters (eo = 1.68 
GeV/fm, ^0 = 1-5 GeV and Tq = 0.2 fm/c) as determined 
from a simultaneous fit to the suppression of both single 
and dihadron spectra [6], except when specified. 

For the study of 7-hadron spectra we follow Ref. [7] 
and define the 7-triggered FF as. 
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where zt = PtIPt — 7r/2 < < 7r/2 is the azimuthal 
angle between the triggered 7 and the associated hadron 
on the away-side. The above 7-triggered FF is a sum of 
FF's of the away-side jets (quark and gluon) weighted 
with the fractional 7-jet production cross sections and 
convoluted with the transverse momentum smearing due 
to NLO processes, which can be exploited to explore dif- 
ferent limit of jet quenching as we will discuss later. The 
hadron-triggered FF's have been similarly defined for the 
study of dihadron spectra [BJ (TB] by replacing photons 
with triggered hadrons in the above equation. 

Shown in the upper panel of Fig. [T] are the calculated 
7-triggered FF's in p + p collisions at the RHIC energy 
as compared to PHENIX preliminary data. They agree 
very well for different values of p^'^. In the lower panel 
of Fig. [1] we show the 7-triggered FF in central Au + Au 
(solid) for 8 < p^'® < 16 GeV/c as compared to p + p 
(dashed) collisions. With the same energy loss parame- 
ters as determined by the single and dihadron suppres- 
sion [6], the NLO pQCD results agree with the STAR 
preliminary data very well. Also shown are the calcu- 
lated hadron-triggered FF's as compared to the experi- 
mental data, which are larger than the 7-triggered FF's. 
This is mainly because the fraction of gluon jets associ- 
ated with a hadron trigger at this range of p^'^ is larger 
than 7-triggered jets and the hadron yield of gluon jets 
is larger than that of quarks. Note that the average ini- 
tial jet energy associated with a trigger hadron is larger 
than that associated with a direct photon for the same 
value of Py'^ due to parton fragmentation. We also show 
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FIG. 1: (color online). 7-triggered and hadron-triggered FF 
's in p + p and central Au-\- Au collisions at the RHIC energy. 
The preliminary data are from [211 122| . Systematic errors of 
experimental data are shown as shaded bars when available. 



in the lower panel the uncertainty, mainly from the scale 
dependence of the FF's, of NLO pQCD results ioi p + p 
due to the choice of factorization scale p,. Here M is the 
invariant mass of the dihadron. We will use p = 0.5p^'^ 
for 7-hadron spectra in this paper unless specified. 

To quantify the suppression of hadron and 7-triggered 
FF 's in central Au + Au relative to p + p collisions due to 
jet quenching, as seen in Fig. [l] one defines the nuclear 
modification or suppression factor. 



Iaa — Daa/Dp 
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for both hadron and 7-triggered FF. Shown in the upper 
panel of Fig. |2] are the calculated nuclear modification 
factors both in LO (dashed) and NLO (solid) for different 
values of the energy loss parameter eq- As compared to 
the same study of the hadron-triggered FF [B], Iaa{zt) 
for 7-triggered hadron spectra at small zt < 0.6, which 
are dominated by volume emission, is more sensitive to 
the energy loss parameter eg and therefore provides a bet- 
ter phenomenological constraint on the medium proper- 
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FIG. 2: (color online), (upper panel) LO and NLO pQCD 
results for the suppression factor for 7-triggered hadrons in 
central Au + Au collisions at RHIC energy with different val- 
ues of energy loss parameter eo- (lower panel) NLO results 
for 7-triggered hadrons are compared to dihadron suppression 
together with STAR data |21j. 
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ties when compared to experimental data. 

In the LO pQCD calculation, transverse momentum of 
the associated jet is balanced exactly by the direct photon 
in tree 2^2 processes. This limits = Pt/Pt — 1- 
NLO, however, the initial jet energy can exceed pj, due to 
radiative correction or broadening in the initial state and 
therefore leads to hadrons with > or z-r > 1. In 
this region, the 7-triggered FF is mainly determined by 
the tail of the radiative broadening which falls sharply 
as a function of the jet transverse momentum. There- 
fore, contributions to the final associated hadron spectra 
Daa{zt) {zt > 1) from 7-triggered jets with even a small 
amount of energy loss will be suppressed. Only those jets 
from surface emission that escape from the medium with- 
out energy loss will contribute, whose FF's are the same 
as in the vacuum. Therefore, the nuclear modification 
factor Iaa{zt) in this region is mainly determined by the 
thickness of the corona of the surface emission. On the 
other hand, jets that have lost finite amount of energy 
before fragmenting into hadrons will contribute to the 
7-triggered FF in the region zt < 0.6 where the nuclear 
modification factor Iaa^zt) is controlled by volume emis- 
sion of jets and is therefore more sensitive to the variation 
of the energy loss parameter. The value and zt depen- 
dence of Iaa{zt) in 0.6 < zt < lA are determined by the 
competition of the two mechanisms of hadron emission. 
One can, therefore, determine the jet energy loss param- 



FIG. 3: (color online). Transverse spatial distributions of the 
initial 7-jet production vertexes that contribute to the final 
observed 7-hadron pairs along a given direction (arrows) with 
Zt ~ 0.9 (upper panel) and zt ~ 0.3 (lower panel). 



eter from different spatial regions in heavy-ion collisions 
by measuring the effective 7-triggered FF in the whole 
range of zt, possibly also for different orientation of the 
7-hadron pair with respect to the reaction plane, achiev- 
ing a true tomographic study of the dense medium. For 
precision studies one should also consider the effect of in- 
trinsic transverse momentum broadening via systematic 
analysis oi p + p and p{d) + A collisions [7 . 

The suppression factors Iaa for both hadron (dashed) 
and 7-triggered (solid) FF's in central Au + Au collisions 
at the RHIC energy are compared with the STAR data in 
the lower panel of Fig.[2j The stronger dependence of Iaa 
for 7-hadrons on Zt is due to the dominance of volume 
emission at small zt- The similarity in value between 
Iaa for dihadron and 7-hadron spectra at z^ « 0.4— 1.0, 
despite of their different trigger biases, is partially due 
to the competition between the larger gluon jet fraction 
(bigger energy loss) and larger initial jet energy (more 
penetration) associated with a hadron trigger. 

To further illustrate the above picture of volume and 
surface emission of 7-triggered jets and their contribu- 
tions to the effective 7-triggered FF at different values of 
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FIG. 4: (color online). A^part dependence of Iaa for 7- 
triggered FF in Au + Au collisions at the RHIC energy. The 
data are from [21] . in NLO pQCD calculation is chosen as 
= 2, 4, 6, 8, 10 GeV/c, respectively (from top to bottom). 



zt, we plot in Fig. |3]the spatial transverse distribution 
of the initial 7-jet production vertexes that contribute to 
the final 7-hadron pairs with given values of zt- The 
color strength represents the 7-hadron yield from the 
fragmentation of the 7-triggered jets after parton energy 
loss. The 7-hadron yields with arbitrary scale are given 
by Eq. The inserted panels are projections of the 
contour plots onto y-axes with (solid) and without en- 
ergy loss (dashed). For zt ~ I- (upper-panel), contribu- 
tions to the final observed 7-hadron pairs indeed come 
mostly from the surface. Contributions from 7-triggered 
jets from the center or the same side of the trigger pho- 
ton are mostly suppressed. However, for small z^ w 0.3 



(lower-panel), contributions to the effective 7-triggered 
FF are seen to come from mostly the whole volume, ex- 
cept near the surface on the side of the trigger photon. 

These pictures of volume and surface emission of 7- 
hadron pairs in heavy-ion collisions will lead to different 
centrality dependence of the suppression factor Iaa{zt) 
for different values of z^- Shown in Fig. [4] are Iaa for 7- 
triggered hadron spectra as functions of the participant 
number iVpart in -Au + Au collisions at the RHIC en- 
ergy for different values of zt as compared to the STAR 
preliminary data. For small values of zt < 0.6, the 7- 
triggered hadron yield is dominated by volume emission 
and therefore the centrality dependence of Iaa is weaker 
than that in the region zt > 1 where surface emission is 
the dominant production mechanism. 

In summary, high pT 7-hadron correlations are stud- 
ied within the NLO pQCD parton model with modified 
parton FF 's due to jet quenching in high energy heavy- 
ion collisions. We demonstrated that volume (surface) 
emission dominates the 7-triggered hadrons spectra at 
small Zt < 0.6 (large zt > 1) due to the underlying jet 
spectra in the NLO pQCD. This will enable one to ex- 
tract jet quenching parameters from different regions of 
the dense medium via measurement of the nuclear modi- 
fication factor of the 7-triggered FF in the whole kinetic 
region, including zt > I, achieving a true tomographic 
study of the dense medium. 
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